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Abstract 

A new analysis and new Molecular Dynamics (MD) simulations of the measurements 
of the mobility Hi of ions in dense supercritical Ar gas are reported. \n shows 
a marked dependence on the distance from the critical temperature T c . A mobility 
defect appears as a function of the gas density and its maximum value occurs below 
the critical density. The locus of points of maximum mobility defect in the P — T 
plane appears on the extrapolation of the coexistence curve into the single-phase 
region. MD simulations quantitatively reproduce the mobility defect near T c . 
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1 Introduction 



Negative ions in supercritical gases are an interesting subject for several rea- 
sons. The formation of anions as a consequence of low-energy electron impact 
to a molecule has important consequences in the chemistry of atmosphere, 
of discharges, and in electro- and biochemistry. Electron attachment and de- 
tachment processes, leading to the formation and decay of anions are strongly 
influenced by the environment pQ. Upon creation and stabilization, they in- 
teract with the environment leading to states that cannot be adiabaticaliy 
obtained by a simple addition of the ion to the surrounding system [2113] . 



The investigation of the transport properties of ions in a dense supercritical 
gas gives important pieces of information about their structure and interaction 
with the host [I] and bridge the gap between dilute gases [5] and liquids [6]. 
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Anions are difficult to produce even in a dense environment, with the notice- 
able exception of oxygen and SF 6 that are very electronegative. In particular, 
2 is an ubiquitous species, present as impurity even in the best purified rare 
gas. Resonant electron attachment |7], followed by stabilization induced by 
three-body collisions j8], produces stable 2 ions that can be easily drifted 
through the gas under the action of an externally applied electric field. 

The O2 ion couples with the gas by means of electrostriction [9]. Its strong 
field polarizes the surrounding atoms which are attracted toward the ion that 
becomes solvated in a cluster. Electrostriction thus induces rapidly varying 
density and viscosity profiles around the ion [4"lllO|lll] that affect the ion mo- 
bility /ij. In these conditions, the ion transport is actually determined by the 
hydrodynamic interaction of the large solvation structure surrounding the ion 
with the gas and is quite insensitive to the ionic species. 

Whereas fii in He at T = 77.2 K [12] does not show any particular features 
as a function of the gas density N up to moderate values, it shows a smooth 
change of behavior when N is varied across the critical density N c in Ne gas 
0.6 K above the critical temperature T c [1]. The dependence of /ij on N for 
N > N c has been described in terms of the Stokes' hydrodynamic formula 
relating ^ to ion radius R and fluid viscosity 77 : \ii — e/6in]R, in which 7] is 
evaluated for the density of the maximum of the density profile induced by 
electrostriction [13J. 

Further measurements in near critical Ar gas [H] have confirmed that the 
coupling of the ion motion with the critical fluctuations via electrostriction 
is responsible for the change of behavior of /ij near the critical point. The 
behavior close to T c has been explained by modifying the Stokes's formula in 
order to account for the fact that the large gas compressibility forces a layer 
of correlated fluid to stick around the ion. In order to test this interpretation 
at microscopic level, Molecular Dynamics (MD) simulations are performed. 



2 Experimental results and MD simulations 

The experimental results of the measurements of the ion density- normalized 
mobility in dense supercritical Ar for T = 157, 154, and 151.5 K are shown 
in fig. [TJ Details on the experiment and on the previous analysis are found in 
literature [TJ]. For a density N m < N c = 8.08 nm -3 \Xi shows a defect whose 
strength increases when T —>■ T c = 150.86 K. The solid line through the data 
on the isotherm closest to T c is obtained by using the Stokes' formula with 
an ion radius R that is enhanced with respect to the radius in absence of 
electrostriction and criticality by an amount that depends on the local value 
of the fluid correlation lenght £ [14J . The good agreement confirms the picture 
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Fig. 1. /UjiV vs N. Open circles: T = 157 K. Squares: T = 154 K. Closed circles: 
T = 151.5 K. Solid line: model prediction. 

of a correlated fluid layer dragged by the ion because of electrostriction. For 
N < 4 nm~ 3 the hydrodynamic description fails to reproduce the data because 
the ion transport leaves the hydrodynamic regime to enter the kinetic one. 

A new analysis of the data reveals that the locus of the maximum mobility 
defect lies on the extrapolation of the coexistence curve into the single phase 
region, as shown in fig. [2J The same occurs for 3 He + ions in 3 He near the 
critical point for T c rs 3.3 K [To] .The gas appears reminiscent of the nearby 
coexistence curve even for T > T c . Droplets of superheated liquid are formed 
in the gas because of critical fluctuations. Electrostriction makes the local 
pressure around the ion rise through the extrapolation of the coexistence line 
into the region of superheated liquid that the strong ion field stabilizes ther- 
modynamically. This situation is also similar to the case of 4 He + ions near 
the melting transition in which the A-line is extrapolated into the region of 
supercooled, electrostrictively-stabilized liquid surrounding the ion [16J. 

To get a microscopic insight of the ionic transport process, classical Molecular 
Dynamics (MD) simulations have been carried out in which the Newtonian 
equations of motion are numerically integrated for suitable interaction poten- 
tials by using freely available Fortran codes [IT]. A system of 100 Ar atoms 
plus 1 ion has been simulated as a function of N for T = 151.5 K. The 
Ar-Ar interaction potential is found in literature [TS|. The dominant long- 
range part of the 07— Ar potential can be calculated by the knowledge of the 
polarizabilities of 07 and Ar and of the quadrupole moment of Ar [19]. The 
repulsive part of the 0.7— Ar potential has been modeled as the interaction of 
the neutral species with an infinite wall endowed with a local charge density 
equal to the electronic charge density of the ion [20] . 
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Fig. 2. Locus of points of minimum ion mobility. Solid line: coexistence curve |21| . 

The equations of motion are integrated by using the Verlet algorithm [I7J 
in time steps of r = lOOrj = 2.42 • 10~ 15 s, where r, = 47ie a h/e 2 is the 
timescale in atomic units (a.u.). After an initial equilibration time, the system 
evolves freely. The coordinates of all particles are recorded every time interval 
t n = (500r)n, n = 1, . . . rif, spanning a total simulation time ~ 2.4 ns for rif = 
2000. In fig. [3] we report a typical example of the mean square displacement 
(r 2 ) of the ion and of Ar for N = 6.5 nm~ 3 and T = 151.5 K. For Ar, it is less 
noisy because of the higher statistical significance due to the averaging over 
100 atoms. As expected, the ion is less mobile than the neutral atoms because 
the enhancement of its interaction with atoms due to its electric field. 

The diffusion coefficient D is obtained as D = lim^t-^oo (r 2 (At)) /QAt. In order 
to improve the statistical significance of the data, a running average of the 
diffusion coefficient _D ra is first calculated by evaluating (r 2 ) for each At by 
moving the time origin across the whole simulation and averaging over the 
results before averaging over the sampling window At. In fig. H]a typical result 
for the running average of _D ra is plotted as a function of the sampling window 
At. The error bars are the results of averaging over the different time origins. 
For practical purposes, the limit in the definition of D is replaced by the 
mean of the running average _D ra . Finally, the Nernst-Einstein relationship 
D/ni = ksT/e is used to calculate the ion mobility. 

In fig. [5] the MD results for T = 151.5 K are compared with the experimental 
data. The quantitative agreement between data and simulation is very good. 
The simulated ^ shows the defect exactly at the same density of the experi- 
ment. The density width of the mobility defect calculated in the simulations 
is larger than the experimental one. It is believed that this effect is due to 
the fact that the temperature during simulation is allowed to substantially 
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Fig. 4. D Ia vs At. 
fluctuate in order for energy to be constant. 

Disappointingly, however, MD simulations do not enlighten the physical rea- 
sons of their success. Actually, MD calculations also provide the ion-atom pair 
distribution function g(r) that gives important information on the structure 
of the ion environment. In particular, the excess number N exc of atoms sur- 
rounding the ion can be calculated. In fig. [H] N exc is reported as a function of 
the distance r from the ion for N close to that of the mobility minimum. The 
formation of the first solvation shell is clearly indicated by the plateau starting 
at r fa 4 A, while a diffuse second solvation shell is observed for r « 9 A. In 
fig. El the number of excess atoms contained in the 1st and 2nd solvation shells 
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Fig. 5. Comparison between experiment (closed points) and MD simulations (open 
points). The error bars are not shown for graphical clarity. Solid lines: eyeguides. 
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Fig. 6. iVexc vs r for N = 6.5 nm~ 3 and T = 151.5 K (left- and bottom axes). iV C xc 
vs N for T = 151.5 K. Squares: 1st solvation shell. Circles: 2nd solvation shell (top- 
and right axes). 



is also shown as a function of N. Interestingly enough, N cxc in the 2nd solva- 
tion shell is strongly peaked for N < iV c , thus suggesting that electrostriction 
efficiently modifies the ion environment when the unperturbed gas density is 
well below iV c . On the other hand, nothing anomalous shows up around the 
density of the maximum mobility defect. We believe that carrying out more 
MD simulations will improve their statistical significance and thus shed light 
on the physical picture beneath the experiment. 
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